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ABSTRACT 

We re-examine the effect of Lyman continuum (A < 912 A) extinction (LCE) by dust in H ii regions 
in detail and discuss how it affects the estimation of the global star formation rate (SFR) of galaxies. 
To clarify the first issue, we establish two independent methods for estimating a parameter of LCE (/), 
which is defined as the fraction of Lyman continuum photons contributing to hydrogen ionization in 
an H II region. One of those methods determines / from the set of Lyman continuum flux, electron 
density and metallicity. In the framework of this method, as the metallicity and/or the Lyman photon 
flux increase, / is found to decrease. The other method determines / from the ratio of infrared flux 
to Lyman continuum flux. Importantly, we show that / < 0.5 via both methods in many H ii regions 
of the Galaxy. Thus, it establishes that dust in such H ii regions absorbs significant amount of Lyman 
continuum photons directly. To examine the second issue, we approximate / to a function of only the 
dust-to-gas mass ratio (i.e., metallicity), assuming a parameter fit for the Galactic H ii regions. We find 
that a characteristic /, which is defined as / averaged over a galaxy-wide scale, is 0.3 for the nearby spiral 
galaxies. This relatively small / indicates that a typical increment factor due to LCE for estimating the 
global SFR (1//) is large (~ 3) for the nearby spiral galaxies. Therefore, we conclude that the effect of 
LCE is not negligible relative to other uncertainties of estimating the SFR of galaxies. 

Subject headings: dust, extinction — H ii regions — galaxies: ISM — infrared: ISM: continuum — 
radio continuum: ISM — stars: formation 
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1. INTRODUCTION 

When we estimate the present star formation rate (SFR) 
in galaxies, we use the luminosity of radiation from young 
massive stars as its indicator. Various observational quan- 
tities are adopted in order to count the number of photons 
originating from these stars. Indeed, we use the data of 
hydrogen recombination lines, ultraviolet (UV), infrared 
(IR), radio, etc., as the indicators of the photon flux and 
the SFR (e.g., Kennicutt 1998). Observational evidence 
suggests, however, that star-forming regions are often as- 
sociated with dust (e.g.. Glass 1999, p. 125). Thus, we 
cannot obtain the true SFR of galaxies unless we correct 
observational data for dust extinction. In other words, 
understanding the extinction property is important in es- 
timating the SFR (e.g., Madau, Pozzetti, & Dickinson 
1998). Therefore, we should examine the properties of 
dust extinction in star-forming regions in detail. 

The extinction by the interstellar dust from UV (A > 
912 A) to near-IR (NIR) has been studied well to date 
(e.g., interstellar extinction curve: Savage & Mathis 1979; 
Seaton 1979; Calzetti, Kinney, & Storchi-Bergmann 1994; 
Gordon et al. 2000). According to those researches, it is 
widely accepted that when we estimate the SFR from the 
Ha luminosity, for example, we should correct decrement 
of the line luminosity due to the UV-NIR extinction. It 
can be performed by using a proper extinction curve and 
the observation of the Balmer decrement (e.g., Osterbrock 
1989). On the other hand, how about the effect of the 
extinction in the Lyman continuum band? Because of the 



decrement of the stellar radiation owing to the Lyman con- 
tinuum extinction (LCE), we may underestimate the SFR 
of galaxies. In this paper, this point is extensively exam- 
ined. 

Indeed, many observations of H ii regions have revealed 
that dust grains are associated with these regions (e.g., 
Ishida & Kawajiri 1968; Harper & Low 1971; Wynn- 
WiUiams & Beckhn 1974; Frey et al. 1979; Mizuno 1982). 
This means a very important fact that the Lyman contin- 
uum photons from the central young massive stars suffer 
extinction by dust grains within H ii regions. As a result, 
the number of Lyman continuum photons contributing to 
hydrogen ionization is reduced. From the intensity of re- 
combination lines or thermal radio radiation, we can only 
estimate the number of the photons really contributing to 
the ionization. Therefore, we inevitably underestimate the 
number of Lyman continuum photons and the SFR unless 
wc correct the observational data for the LCE by dust. 

In studies on each individual H ii region, it has been 
reported that a part of Lyman continuum photons are ab- 
sorbed by dust if H ii regions contain dust grains (e.g., Pet- 
rosian. Silk, & Field 1972; Panagia 1974; Natta & Panagia 
1976; Sarazin 1977; Spitzer 1978; Mathis 1986; Aannestad 
1989; Shields & Kennicutt 1995; Bottorff et al. 1998). Ac- 
cording to Petrosian et al. (1972), for instance, the fraction 
of Lyman continuum photons used by hydrogen ionization 
is estimated to be 0.26 for the Orion nebula. Their result 
has shown that the chance which the Lyman continuum 
photons are extinguished by dust is significantly large. On 
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a galactic scale, Smith, Biermann, & Mezger (1978) have 
examined LCE by dust for a large number of radio H ii 
regions in the Galaxy, and then, have shown about half 
Lyman continuum photons from exciting stars in many 
H II regions arc absorbed by dust within these regions. 
However, they did not discuss the effect of the LCE on 
estimating the SFR of galaxies. 

Unfortunately, it seems that we often estimate the SFR 
of galaxies from Lyman continuum flux (e.g., Kennicutt 
1983) without any correction of the LCE effect. Hence, it 
is indispensable for us to re-examine LCE by dust within 
H II regions and to discuss the effect of LCE on deter- 
mining the SFR of galaxies. First, we formulate the way 
of estimating the fraction of Lyman continuum photons 
used by hydrogen ionization, and determine the fraction 
in each individual H II region of the Galaxy in section 2. 
We also present another way for estimating LCE by dust 
in section 3. Then, we discuss the effect of LCE by dust 
on estimating the galaxy-scale SFR in section 4. Finally, 
our conclusions are summarized in the last section. 

2. FRACTION OF LYMAN CONTINUUM PHOTONS 

CONTRIBUTING TO HYDROGEN IONIZATION 

The Lyman continuum photons from young massive 
stars in an H II region are absorbed by dust grains as 
well as by neutral hydrogen atoms. Then, the number 
of Lyman continuum photons contributing to hydrogen 
ionization, -/V^y, is smaller than the intrinsic number of 
Lyman continuum photons, Nlj. We parameterize this 
effect by /, which is defined as the fraction of Lyman con- 
tinuum photons contributing to the hydrogen ionization; 
iV^y = /iVLy. Unfortunately, we are able to estimate only 
N^y from observations of recombination lines or thermal 
radio continuum. This means that the "real" SFR, which 
should be estimated from A^Ly, is 1// times larger than 
the "apparent" SFR estimated from N^^y Hence, if the 
parameter, /, is much smaller than unity, the effect of 
LCE by dust becomes very important to the estimation of 
the SFR. We call 1// correction factor or increment fac- 
tor for the SFR. In this section, we examine the value of 
/ in each individual H II region in the Galaxy to find the 
correction factor quantitatively. 

2.1. Formulation 

First, when no dust grains exist in an H II region, the ra- 
dius of the region is estimated to be the Stromgren radius, 
rg. According to Spitzer (1978), we express the number 
of Lyman continuum photons emitted per unit time, A^Lyj 
by using the Stromgren radius, rs- 

Nlj = y rs^nenpa(2) , (1) 

where rie and rip are the number densities of electrons and 
protons, respectively, and a^^^ is the recombination coef- 
ficient excluding captures to the n = 1 level. Here, we 
employ the Case B approximation, which is the assump- 
tion that H II regions are optically thick for the photons of 
all Lyman-emission lines (e.g., Osterbrock 1989). We also 
assume spherical symmetry and spatial uniformity of H II 
regions in this paper for simplicity. 

Next, we consider the effect of LCE by dust on the size 
of H II regions. The actual radius of the ionized region be- 
comes smaller than the Stromgren radius in equation (1) 



by the LCE effect. When we express the actual ionized 
radius as ri, we obtain 



fNLy 



3 ('21 



(2) 



Moreover, if the ratio of n to rs is expressed by j/i, we 
obtain the following relation from equations (1) and (2), 



f = V?- 



(3) 



To estimate /, we need to find the quantitative relation 
between / and the amount of dust. We determine / in 
terms of the optical depth of dust for Lyman continuum 
photons. According to Hirashita ct al. (2001), then, we 
define Ts,d as the optical depth of dust for Lyman con- 
tinuum photons over the Stromgren radius in order to 
formulate the dependence of / on the dust-to-gas ratio. 
Here, rs,d is approximated to be the optical depth at the 
Lyman limit (912 A). When we adopt the Galactic extinc- 
tion curve, the dust extinction at 912 A is about ViEs-v 
mag. Thus, Ts,d — 13i?B_v/2.5 loge ~ V1Eb-v- If we 
assume that the dust-to-gas mass ratio, I?, is proportional 
to Eb-v /N-n, where A'h denotes the column number den- 
sity of hydrogen, then, according to Spitzer (1978) and 
Hirashita et al. (2001), 
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where V is estimated with respect to a typical Galactic 
value, 6 x IQ-^. Thus, 
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Here, we define A^h as the column density over the 
Stromgren radius. That is, A^h = nnrs, where nn de- 
notes the volume number density of hydrogen. Therefore, 
equation (5) is reduced to 
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(6) 

where rs is eliminated by equation (1), and we suppose 
that Up !v He !v riH (fully ionized in the ionized region) 
and that the temperature of H ii regions is 10'' K. Since 
wc cannot observe the real number of Lyman continuum 
photons, A^Ly, we rewrite equation (6) in terms of the ap- 
parent one, Nly{= yi^Ni^y); 
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(7) 

In addition, Spitzer (1978) derived the following relation 
between and rs^d^ 



3 / y'^e^^^-'^dy = 1 
Jo 



(8) 



where y is the radius normalized by rs- Integrating equa- 
tion (8), we obtain 



rs,d' = 3{e^<»(Ti-2Td + 2)-2}, 



(9) 



where rd is defined as yirs,d- It is worthwhile to note 
that Td is identified with the optical depth of dust over 
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the actual ionized radius, n. According to equation (3), 
f = V"? = (7'd/7's,d)^- Thus, we finally obtain 



/ = 



3{e-<i(r2-2Td + 2)-2}" 



(10) 



This is the same as equation (8) in Petrosian et al. (1972). 

Now, once nn (~ ne)^ N^^, and V are determined, we 
derive Td(= yirs.d) via equation (7), and then, we cal- 
culate / by using equation (10). In Figure 1, we show 
the / vs. Td relation obtained from equation (10). Which 
quantity determines / the most effectively? According to 
figure 1, decrement of / is due to increment of the dust 
optical depth. Thus, to find the most effective parame- 
ter for /, let us examine the parameter dependence of the 
optical depth. We find in equation (7) that the optical 
depth increases according to the dust-to-gas ratio, Lyman 
continuum flux, and density of Hii regions. Then, when 
any of those three parametes increase, / decreases. By 
the way, among the power indices in equation (7), that of 
the dust-to-gas ratio is the largest. Then, we find that / 
depends most largely on the dust-to-gas ratio. 

One might think that the effect of helium on the ion- 
ization structure of H II regions is important (e.g., Os- 
terbrock 1989). In H ii regions, there arc extremely en- 
ergetic photons whose wavelengths are shorter than 504 
A. Such photons can ionize neutral helium. Hence, the 
number of Lyman continuum photons absorbed by hydro- 
gen seems to decrease because of the absorption by helium. 
However, fortunately, almost all (96% according to Mathis 
1971) photons produced by helium recombination can ion- 
ize neutral hydrogen. Therefore, we can assume safely that 
the number of Lyman continuum photons does not signifi- 
cantly decrease owing to the presence of helium. Moreover, 
the fraction of the photons with shorter wavelengths than 
504 A is small (0.14 in our calculation adopting the stellar 
spectra assumed to be the Plank function and Salpeter's 
IMF). Thus, the effect of helium on the photon count is 
small in estimating the SFR. Another effect of helium ap- 
pears on the electron number density. For example, in the 
derivation of equation (6), we must replace Uq = nu with 
Tie = nu + nue- Since helium abundance is typically 0.1 
by number, Uc = Linn- Hence, we consider this effect to 
be less significant than that caused by some other uncer- 
tainties, e.g., the change of the IMF. Therefore, we neglect 
these effects of helium in this paper. 

2.2. Value of f for Individual H ii Region in the Galaxy 

Let us estimate / for the Galactic H II regions in order 
to estimate the effect of LCE by dust on determining the 
SFR. To do this, as formulated in the previous subsection, 
we need the data set of the number of Lyman continuum 
photons, the electron (or hydrogen) number density, and 
the dust-to-gas ratio. 

The amount of dust is considered to be related to the 
metallicity. In this section, we determine the dust-to-gas 
mass ratio, V, from the observational metallicity (abun- 
dance of oxygen) for each H II region, using a global re- 
lation between T> and (0/H) proposed and modeled by 
Hirashita (1999a, b).'^ We display in figure 2 the global 

^ The parameters are set as /in,o =0.1 and /3aoc = 2/3sN = 10. Here, /in 
/3acc and /3sn are defined as gas consumption timescale (gas mass divided 
dust destruction timescale, respectively (see Hirashita 1999a for details). 



relation between V and 12 -I- log(0/H). This is the same 
as Figure 4 of Hirashita et al. (2001). Here, we adopt 
12 + log(0/H) = 8.93 as the solar abundance (Cox 2000). 
We should note here that this 'D-(0/H) relation is the av- 
erage relation for the interstellar medium (ISM) in nearby 
spirals and dwarfs. That is, we assume that V-{0/B.) re- 
lation of H II regions is not different significantly from the 
averaged I?-(0/H) of ISM. This assumption may not be 
very good, but the dust properties in H II regions are still 
uncertain. Hence, as a first step, we apply the I?— (0/H) 
relation of the ISM to each H ii region. 

We make two sets of the suitable samples to estimate 
/ of each individual H ii region. One is the sample of 
seven representative giant H ii regions in the Galaxy. Their 
properties are summarized in Table 1. These H II regions 
are famous and have been studied well. Thus, the uncer- 
tainties of their data, especially their distance (i.e., their 
luminosity), can be regarded as small. Their apparent 
number flux of Lyman continuum photons (-^Ly) are esti- 
mated from radio observations and more than about 10^® 
s""'^. Hereafter, we call this sample luminous sample. The 
determined / values are also shown in column (7) of Ta- 
ble 1. We find a good agreement between Petrosian et al.'s 
/ and ours for the Orion nebula (/ = 0.26). The mean / 
of this sample is 0.35. 

The other sample is a part of H ii regions observed and 
investigated by Caplan et al. (2000) and Deharveng et al. 
(2000). We select the H ii regions whose kinematic dis- 
tance and photometric distance are consistent within a 
factor of 1.5 from their sample. This is because the un- 
certainty in distance affects the estimation of the number 
of Lyman continuum photons and then prevents us from 
determining / reasonably. These H ii regions tend to be 
less luminous than the sample of Table 1 and their ap- 
parent number flux of Lyman continuum photons (-^Ly) 
are estimated from their Ha luminosities. Hereafter, the 
sample is called to less luminous sample. Their proper- 
ties are tabulated in Table 2. Caplan et al. (2000) also 
observed the Orion nebula and M16, both of which are 
also included in the sample in Table 1, but their Ho; lu- 
minosities are smaller than that expected from the radio 
observations adopted in Table 1. This is because Caplan 
et al. (2000) did not observe whole parts of these nebu- 
lae. Thus, we exclude the Orion nebula and M16 from the 
sample of Table 2. The finaly determined / of the sample 
distributes over 0.5-1, and the mean is 0.8. 

In Figures 3 (a)-(f), we present several relations among 
various parameters of our sample H ii regions in Tables 1 
and 2. Opcni symbols are H II regions in Table 1, and 
filled symbols are those in Table 2. In the panels of (a)- 
(c), we show some relations among the basic parameters to 
determine /. From these panels, we can check our stand- 
point for the subsequent considerations. Firstly, we may 
find, in panel (a), the sample regions with higher metal- 
licity have larger observed Lyman continuum flux, A^Ly 
Indeed, the sample correlation coefficient is 0.41, so that 
the correlation exists with a confidence level of 95% for 
the sample. However, we find no correlation for the sam- 
ple of Table 2 alone. We also cannot find the correlation 

, O is the dust mass fraction in the material injected from stars, and 
by star formation rate) normalized by dust growth timescale and by 
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for the extragalactic H ii regions.'' Thus, we consider that 
the correlation in panel (a) may be caused by our sam- 
ple selection. Hence, we assume it conservatively in this 
paper that iV^y dose not depend on the metallicity. In 
panel (b), the correlation between the metallicity and the 
electron number density, rie, is not found. Thus, we dose 
not consider that rig depends on the metallicity. In panel 

(c) , interestingly, wc can find a correlation between and 
N^y for only samples in Table 2. 

Next, in the panels of (d)-(f), we present the dependece 
of / on the parameters of H ii regions. The estimated / is 
regarded as a function of the metallicity in panel (d), i.e., 
as the metallicity in an H ii region increases, / becomes 
smaller. This is because, in our formulation, the dust- to- 
gas ratio (i.e., metallicity) has a positive dependence on 
the optical depth of dust as shown in equation (7), and 
then, a larger dust content leads to a larger optical depth 
of dust and smaller /. Although the correlation in (d) is 
not the result directly derived via the analysis of the ob- 
servational quantities, we confirm a trend that / depends 
on the mctalicity (see also Hirashita et al. 2001). In panel 
(e), we also find / is related to N^^^ which has a positive 
dependence on the dust optical depth (eq. [7]). This issue 
will be discussed again in section 3. We cannot find a sig- 
nificant correlation in panel (f), probably because ric has 
a weak dependence in the dust optical depth (index 1/3 
in eq. [7]) and its dynamic range is small relative to N[^y 
which has the same dependence in equation (7). 

Now, we concentrate on the dependence of the dust-to- 
gas ratio (i.e., metallicity) on /, because the variation of 
dust-to-gas ratio has the largest influence on determining 
the dust optical depth as seen in equation (7). That is, 
we assume that the optical depth of dust is a function 
of the dust-to-gas ratio only. If we consider AT^^ and rie 
dose not depend on V (or metallicity), we can re-express 
equation (7) as 

Td ^ yirs,d = X (^^^) , (11) 

where ra is the dust optical depth over the actual ionized 
radius (see section 2.1), and x is a factor containing the 
dependence of iV^y and nH(~ ric). We find that the best 
fit value of X for the adopted Galactic H ii regions is 1.8, 
although the dispersion is large. In Figure 4, we display 
the estimated / of the individual H ii region in our sample 
as a function of its dust-to-gas mass ratio, V. The open 
and filled symbols represent the sample of Tables 1 and 2, 
respectively. This figure is basically identical to Figure 3 

(d) . The sohd line is the best fit line {x — 1.8). We also 
show the lines of a; = 1.0, and 3.0 for comparisons. 

Here, let us discuss the increment factor for the SFR, 
1//, which is introduced at the beginning of section 2. It 
is useful for us and readers to check quantitatively how / 
is affected by "D. As a representative case, we consider the 
model of x = 2, which corresponds with A'^y ~ lO^'^s"^ 

and rie — 100 cm""^ (nearly equal to the Lyman continuum 
number flux and electron number density of the Orion neb- 
ula). Numerically, from the set of equations (7) and (10), 
we find / to be about 0.2 for D = 6 x 10~^ (the typical 
value of the ISM in the Galaxy; Spitzer 1978). If X> is 1/3, 



1/10, or 1/100 of the Galactic value, / becomes 0.6, 0.9, 
or nearly unity, respectively. Thus, it is found that the 
increment factor for the SFR, 1//, increases from almost 
unity to about 5 as increases from 6 x 10"** (1/10 of a 
typical Galactic value) to 6 x 10""^. Clearly, the correction 
factor for the SFR is very sensitive to the dust-to-gas ratio 
(or metallicity), especially around T> ~ 10~^. Hence, we 
should determine X> precisely to find an actual / and SFR. 
In any case, we must take account of the effect of LGE by 
dust at least for the objects whose metallicity is as high 
as that of the Orion nebula. 

Smith et al. (1978) have also determined / for a nu- 
merous radio sample of giant H ii regions in the Galaxy. 
Although our equation (7) is equivalent to their equa- 
tion (A. 2), there is a difference in the method of deter- 
mining the dust-to-gas ratio of H ii regions. Indeed, they 
used an empirical relation of the absorption cross section 
as a function of the Galactocentric radius (Churchwell et 
al. 1978), whereas we determine T) for each H ii region 
individually as described above. The determined mean / 
for the giant H II regions by Smith et al. (1978) is 0.56. 
Since N'-^^ of their sample is larger than 10^^ s~^, their 
sampk^ H II regions correspond to our sample in Table 1, 
which has mean / = 0.35. The estimated / values by us 
may be systematically small. This is due to the differences 
of some adopted parameters in calculation, for example, 
Ntil Eb^v ■ Also, Smith et al. (1978) took account of the 
effect of helium which we neglect as mentioned in the pre- 
vious subsection. The parameter of / in Smith et al. is 
thus determined as the sum of the fractions of Lyman con- 
tinuum photons absorbed by both hydrogen and helium. 
This increment by helium in / is about several % (up to 
10%) in their framework. If the increment is removed. 
Smith's mean / will approach our value. 

Throughout this section, we have found that the deter- 
mination of T) is very important to estimate the "real" SFR 
(or /) from an observational count of Lyman continuum 
photons. To obtain D accurately, it is a recommendable 
way that we use the luminosity of IR thermal emission of 
dust in H II regions. Since dust grains absorb the radiative 
energy from the young massive stars and then re-emit the 
energy in the IR range, the amount of dust is determined 
from the IR luminosity. In the next section, we discuss 
LCE by dust from observation of IR luminosity. 



3. CORRELATION BETWEEN IR AND Ha/RADIO 
LUMINOSITIES FOR H II REGIONS 

In the previous section, we estimate the amount of dust 
(the dust-to-gas ratio) in H ii regions from their metallic- 
ity. Since the amount of dust is also reflected in the IR 
emission, / is also determined from IR luminosity inde- 
pendently. In this section, then, we determine / by using 
the IR luminosity of H II regions in order to check consis- 
tency between the results of the previous and the current 
section. 

Based on the theory of the IR emission from dust in 
H II regions by Petrosian et al. (1972), we formulate the 
relation between the observed IR luminosity, Lf^^, and 
the observed number of Lyman continuum photons, A/'^y 

* We examine whether H ii regions with higher metallicity have larger luminosity for the data of H ii regions in MlOl presented by Scowen, 
Dufour, & Hester (1992). The data can be taJsen from ADC, NASA/Goddard or CDS, Strasbourg, Prance via on-line. 
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(Inoue, Hirashita, & Kamaya 2000, 2001). The derived 
formula is 

Lm\S, - 1000/xm)/LQ _ 0.44 - 0.28/ + 0.56e 

7V[y/5.63 X 1043 s-i 7 ' ^ ' 

where e is the average efficiency of the dust absorption 
for UV-optical photons from OB stars. Here, we have 

also adopted the Salpeter's IMF (the stellar mass range 
is 0.1 100 Mq) and proper stellar properties. Moreover, 
the intrinsic number of Lyman continuum photons, Ni^y, 
is 5.63 X 1043 s~^ per unit solar luminosity of the bolomet- 
ric luminosity of OB stars, in our calculation. As defined 
in the previous section, N^^ = fNi^y. Instead of deter- 
mining the SED of the IR radiation, we consider the total 
IR luminosity of dust, Lf^^^, in the range of 8 - 1000 iim 
which covers almost the whole wavelength range of dust 
emission. 

First, let us consider the correlation between the lumi- 
nosities of IR and Ha for H ii regions. We make a cross- 
reference between the sample of Caplan et al. (2000) and 
the IRAS Point Source Catalog {IRAS PSC; Joint IRAS 
Science Working Group 1985). We select the sample H 
II regions by the following criteria; (a) the separation be- 
tween the positions of the Ha observations and sources in 
IRAS PSC is within 100". (b) log(F25/Fi2) > 0.4. (c) 
log(i^6o/^^25) > 0.25. (d) i^ioo > 80 Jy. Here, F12, F25, 
Feo, and Fioo are the lists' fluxes at 12, 25, 60, and 100 
^m, respectively. The criterion (a) is based on the typi- 
cal radius of optical H 11 regions (~ 1-10 pc) and IRAS 
positional error (< 20", Beichman 1985). The criteria (b) 
- (d) are proposed by Hughes & Maclcod (1989). The 
confidence level of the set of criteria for true association 
between IRAS sources and H 11 regions is more than 80% 
(Hughes & Macleod 1989; Codella et al. 1994). As a result, 
we obtain 18 H 11 regions (10 samples are also in Table 2) 
and their properties arc found in Table 3. 

When we suppose that the dust emission can be fit- 
ted by the modified black-body radiation of 30 K (spec- 
tral index is 1), IRAS 40 - 120 /um luminosity, LxraSi 
is L/KA5(40 - 120 Aim) = 0.6L^^"'(8 - 1000 jim). This 
correction factor is consistent with that of Calzetti et al. 
(2000). On the other hand, A^l^ is estimated in terms 
of the Ha luminosity corrected for interstellar extinction: 
Af[Js-i = 2.83 X lO45LHa/i0, where wc assume the 
Case B and the electron temperature of 104 K (Osterbrock 
1989). Thus, equation (12) is reduced to 

LiRAs _ FiRAs _3-^ / 0.44-0.28/ + 0.56e \ ^^^^ 

where Firas and Fua are fluxes of IRAS and Ha, respec- 
tively. This is a theoretical relation between Fjras and 
-Fhq,. Since a flux ratio does not depend on the distance 
to the objects, equation (13) is free from the imcertainty 
in the distance. 

The average efficiency of dust absorption for UV-optical 
photons from OB stars, e, is defined as e = 1 — 1Q~^'^''^^\ 
where {A\) is an average dust absorption over the wave- 
length range of A = 1000 - 4000 A in units of magnitude. 
This range is suitable for our purpose because the number 
of photons of A > 4000 A emitted by OB stars are negligi- 
ble. If we adopt the Galactic extinction curve (Savage & 

^ If the correlation coefficient is more than 0.45 for 18 samples, we can 



Mathis 1979), then we have {Ax) ~ 7.2E{B - V). Thus, 
we obtain the following equation: 

e = 1 - l{)-^E(B-v) (^4) 

Using this equation, wc find e ^ 1 for the H 11 regions in 
Table 3. 

Once we obtain the value of e for each H il region from 
equation (14), we can determine each / from IR and Ha 
fiuxcs by using equation (13). We show the estimated / 
in column (7) of Table 3. By definition, / < 1. However, 
for some sample H II regions, / appears to be greater than 
unity. The reason for the discrepancy will be discussed in 
the last paragraph of this section. 

By the way, we compare the theoretical FjRAs^Fna re- 
lation with observational data on the Fjras vs. Fna plane 
in Figure 5. Adopting e = 1, we make three lines for three 
values of /: the solid, dotted, and dashed lines for / = 1, 
0.5, and 0.1, rc^spcctively. Moreover, we make the dash- 
dotted line for the set of / = 1, e = 0. This corresponds 
to the case that only the Lyman a photons are absorbed 
by dust grains. As long as all Lyman a photons are ab- 
sorbed by dust within H II regions and re-emitted in the 
IR range (Spitzer 1978; Hirashita et al. 2001), the area 
below the dash-dotted line on the Firas vs. Fhq plane 
is the forbbiden area within the framework of our the- 
ory. That is, the dash-dotted line stands for the lower 
boimdary of IRAS flux corresponding to Ha (Lyman con- 
tinuum photon) flux in our model. While the correlation 
between Fjras and is expected theoretically from 
equation (13), we do not find a significant correlation in 
observational data points in figure 5 (the sample correla- 
tion coefficient is 0.35).'' This point is also examined in 
the final part of this section. 

On the other hand, it is well known that there is a tight 
obscrvaticmal correlation between the IR himinosity and 
the number of Lyman continuum photons estimated from 
the radio luminosity (Wynn-Wilhanis & Beckhn 1974). 
Using equation (12), we compare this observational cor- 
relation with the theoretical one in figure 6. Here, we 
read the data directly from the figure in Wynn- Williams 
& Becklin (1974) and re-plotted them. The sample cor- 
relation coefficient of their data is 0.96, showing a stong 
correlation. The parameters of lines are the same as those 
in figure 5. Comparing the data and the model lines, we 
consider the IR luminosity of Wynn- Williams & Becklin 
(1974) (40 - 350 /zm) is nearly equal to the whole dust-lR 
luminosity in equation (12) because the dust luminosity of 
the wavelength range of A < 40 fxm or A > 350 /xm is negli- 
gible (Calzetti et al. 2000). From figure 6, we find that the 
theoretical lines agree with the data of Wynn- Williams & 
Becklin (1974) excellently. This shows that our model is 
supported by the observation. 

If we assume e = 1 for the sample of Wynn- Williams 
& Becklin (1974), we can determine their / from equa- 
tion (12). We find / = 0.45 as an average in the current 
sample. Now, we compare it with another mean / value 
estimated by the method of the previous section. Since H 
II regions in figure 6 have N^^^ > 104^ s"-'^, we shall choose 
18 regions whose Nl^y is greater than 104^ s~^ from Ta- 
bles 1 and 2. Thus, we obtain a mean / as being 0.56 
by the method of section 2. Fortunately, we find a good 
agreement between the two /'s. 

conclude that a correlation exists with a confidence level of 95%. 
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By the way, / values estimated by the both methods 
are not fully consistent with each other. For example, the 
sample of Wynn- Williams & Becklin (1974) includes the 
Orion nebula, M8, and Ml 7 as the sample of Table 1 does. 
Their / values via equation (12) are 0.44, 0.59, and 0.59, 
respectively. Obviously, these values are larger than those 
of the Orion nebula, M8, and Ml 7 in Table 1. This may 
be caused by the aparture difference between two pho- 
tometories, IR and radio (see also the last paragraph in 
this section). This discrepancy may also result from our 
assumption in equation (12) that all of the energy of pho- 
tons not used by hydrogen ionization is converted into IR 
radiation. Generally, the efficiency of conversion to IR 
radiation is never unity because of the escape of Lyman 
continuum photons from H II regions or the scattering of 
these photons by dust. As a result, we expect rather large 
IR luminosities for the H II regions by using equation (12). 

Moreover, we cannot find in Figure 6 the trend that 
the more luminous H ii regions have the smaller /. This 
trend was seen in section 2.1, where we determined / from 
the metallicity. This may be the effect of IR cirrus. The 
observed IR flux originates from dust associated with the 
star forming regions and from diffuse dust heated by the 
interstellar radiation field. If we consider that the IR sur- 
face brightness of the cirrus component is almost constant 
in a galaxy, the effect of cirrus contamination in the ob- 
served IR flux of an H ii region is larger as its IR flux is 
lower. In the case, we may overestimate the IR luminosity 
for low luminosity regions, so that we may underestimate 
/ of such regions. 

Figure 6 has other interesting implications: First, the 
dash-dotted line (f — 1 and e = 0) means that only Ly- 
man a photons heat dust grains and only this energy is 
re-emitted in IR. Obviously the predicted IR luminosity 
in this case falls short of the observed one. Next, the solid 
line (/ = 1 and e = 1) represents the case that only the 
energy of photons whose wavelength is longer than the Ly- 
man limit is absorbed by dust and re-emitted in IR. Even 
in this case, the predicted IR luminosity is not enough to 
explain the observed one. Hence, we conclude that a part 
of Lyman continuum photons ought to be aborbed and 
reprocessed by dust. The fraction of Lyman continuum 
photons absorbed by dust is about 0.5 or more. In terms 
of the fraction of the photons used by hydrogen ionization, 
/ < 0.5 in many H II regions of our Galaxy. This is one of 
the main conclusions in the current paper. 

To close this section, we must discuss a few remain- 
ing questions in the previous paragraphs: Why is / de- 
termined by our theory greater than unity for some H II 
regions in Table 3? Why cannot we find the observational 
correlation between Fjhas and -Fhq in figure 5? These 
may be because we use the IRAS PSC. Since the Galactic 
H II regions lie close to us and extend over large areas, they 
can never be regarded as point-like sources. Hence, the IR 
fiux for H II regions in IRAS PSC may be underestimated 
since the aperture effect of H II regions is not neglected. 
Therefore, IRAS PSC may not be suitable when we mea- 
sure the IR flux of Galactic H ii regions (see also Chan & 
Fich 1995). Conversely, the correlation between the IR lu- 
minosity and the number of Lyman continuum photons in- 
ferred from radio observations is tight. This observational 
correlation is explained very well by our theory. This is be- 



cause the IR luminosities of sample in Wynn- Williams & 
Becklin (1974) are observed in detail individually. Seeing 
the observational tight correlation in figure 6, the absence 
of the correlation in figure 5 is caused by flux underestima- 
tion of IRAS PSC due to the finite aperture effect. Thus, if 
we use the image data of the H li regions observed by IRAS 
or recent/future facilities such as ISO, SOFIA, SIRTF, or 
ASTRO-F (see Takeuchi et al. 1999 for a summary of the 
facilities), the observational correlation between Fir and 
Fua will appear. 

4. EFFECT OF / ON GALACTIC SCALE 

In this section, we estimate / averaged on galactic scale 
for the nearby spiral galaxies. Then we discuss the effect 
of LCE by dust on determining their galaxy-wide SFR. 
In the privious sections, we have discussed LCE by dust 
in H II regions. As shown in section 2, the parameter / 
strongly depends on the dust-to-gas mass ratio. For the 
Galactic H ii regions, we show that x '--^ 2 is appropriate 
in equation (11). That is, the Orion nebula is likely to 
be a typical Galactic H ii region in gas density and Ly- 
man continuum flux. Here, we consider the Orion nebula 
is also a typical H ii region in the other galaxies. In short, 
we assume x = 2 in equation (11) to be universal for the 
galaxies. Then, we estimate a typical / for the nearby 
spiral galaxies from their dust-to-gas ratio. By the way, 
we can determine / of the galaxies from equation (13) as 
perfomed for the Galactic H ii regions in section 3. This 
alternative way is discussed in the last part of this section. 

In order to determine their dust-to-gas ratio, we esti- 
mate the amount of dust mass in the nearby spiral galaxies 
from their IR emission re-processed by dust. As considered 
in the final paragraph of section 3, we do not use IRAS 
PSC since we want to estimate the entire IR luminosity 
of the sample galaxies safely. We adopt only the sample 
galaxies whose entire images by IRAS are available. We 
note another issue here. If we convert the IR emission cal- 
culated only from IRAS data to the dust mass, the total 
amount of dust will be underestimated. This is because 
there is large amount of cold dust (10-20 K) in the spiral 
galaxies, even in those with star burst regions (Calzetti et 
al. 2000 and references therein). Such cold dust cannot be 
traced by IRAS bands. On the other hand, ISO has better 
detectability for the cold dust than IRAS hecause ISO cov- 
ers longer wavelengths than IRAS. Thus, we expect that 
ISO data improve the estimation of dust mass. 

Alton et al. (1998) present the imaging properties of 
seven nearby spiral galaxies obtained by both IRAS and 
ISO. Their sample galaxies are NGC 134, NGC 628 (M74), 
NGC 660, NGC 5194 (M51), NGC 5236 (M83), NGC 6946, 
and NGC 7331. Some basic properties of these galaxies 
are tabulated in tables 1, 2 and 3 of Alton et al. (1998). 
We estimate the dust-to-gas mass ratio for Alton's sample 
galaxies, and then, determine their / averaged over galac- 
tic scale. To avoid the confusion from / of the H ii regions, 
in this section, we express the galactic / as being /. Of 
course, each sample galaxy has unique /. 

We, first, calculate the dust temperature of the sample 
galaxies from the ratio of flux densities at 100 and 200 
fim. In this calculation, we assume the spectrum of dust 
emission to be the modified black-body radiation with the 
spectral index of 1. Here, we note that the determination 
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of the dust temperature strongly affects the successive esti- 
mation of the dust mass. The obtained dust temperatures 
are shown in column (5) of Table 4. Next, we estimate the 
dust mass, Mj, in the sample galaxies from the observed 
flux density at 100 nm via the following equation: 



(15) 



where a, p, and D are the grain radius, grain density, and 
galaxy distance, respectively. F^,, Q^, and BuiT) are the 
observed flux density, grain emissivity, and the Plank func- 
tion of the temperature T at the frequency u, respectively. 
According to Hildebrand (1983), we set (4/3)a/9/(5^ ~ 0.04 
g cm~^ at the frequency of 100 /xm (a = 0.1 /im and p = 3 
g cm~^). The estimated dust masses are tabulated in col- 
umn (6) of Table 4. Then, we determine the dust-to-gas 
mass ratio, 2?, of the sample galaxies by using their to- 
tal (H I + H2) gas masses reported in literatures (column 
7). The values of 2? are in column (8), and the average is 
about 5 X 10~^. This is quite consistent with V estimated 
by Alton et al. (1998). Stickel et al. (2000) also claim that 
the median value of the dust-to-gas mass ratios of sample 
galaxies of ISOPHOT 170 /im Serendipity Survey (Bogim 
et al. 1996) is about 4 x 10~^. Thus, our estimation of V 
is reasonable. 

It still remains open whether 2? in H II regions can be 
approximated to an averaged V of the ISM in the galax- 
ies. As seen in subscx;tion 2.2. there is a large variation of 
V (i.e. observed metallicity) of the Galactic H 11 regions. 
Moreover, most of the sample H 11 regions in Tables 1 and 
2 have smaller V than the typical value for the Galactic 
ISM (6 X 10-3; Spitzcr 1978). It may suggest that X> in H 
II regions is systematically smaller than that of the ISM. 
However, it is reported that the mean V in the ISM is not 
so different from that in H II gas. Sodroski et al. (1997) 
have determined the dust-to-gas mass ratio of the Hi, H2, 
and low-density (no ^ 10 cm^'') H II gases in the Galaxy 
from the data of COBE. According to them, the mean 
value of V in the Galactic ISM (H i -h H2 -h H 11) is 4.6- 
6.6x10-3, and that in H 11 gas is 5.1-10x10-^. With the 
uncertainty of a factor of 2, therefore, we assume that T) 
in H II regions is approximated to the mean V of the ISM 
in the galaxies. 

Here, let us calculate a characteristic value of the /- 
parameters for the sample seven spiral galaxies, adopting 
their dust-to-gas mass ratio (column 8 in Table 4), equa- 
tions (10) and (11) with x = 2. The derived values of / 
are tabulated in column (9) of Table 4. The averaged / is 
0.3. This is nearly equal to the mean / of the sample H 11 
regions in Table 1 which are more luminous than the Orion 
nebula. This small / means that the correction factor of 
the SFR (1//; see also section 2) for these spiral galaxies 
is large (typically ~ 3). 

Although / values of H 11 regions distribute over a wide 
range (0.02-1) as shown in Tables 1 and 2, / is approxi- 
mately equal to the averaged value of the luminous H 11 
regions (i.e., the sample in Table 1). This may be caused 
by the effect of luminosity function of H II regions. Indeed, 
it is implied that the total luminosity of H 11 regions is 

® We assumed implicitly in the previous sections that the IR emission from the H II regions is not contaminated significantly by the cirrus 
emission. This is because the H 11 regions are the star forming regions themselves. However, we may not be able to neglect the contamination 
of cirrus for H II regions with low luminosity. 



dominated by the regions more luminous than the Orion 
nebula (Kennicutt, Edgar, & Hodge 1989; Walterbos & 
Braun 1992; Wyder, Hodge, & Skeffon 1997). Thus, / 
may be shifted toward smaller values because we expect 
smaller / in more luminous H II regions (see section 2.1). 
Therefore, we think it reasonable that a characteristic / 
for nearby spiral galaxies is about 0.3. 

Now we discuss the importance of the increment of a 
factor of 1// ~ 3 for the estimated SFR. There are two 
other major uncertainties for the estimation of the SFR. 
One is related to the initial mass function (IMF). The 
choices of a specific IMF and its upper/lower limit masses 
affect the coefficient of the conversion formula of the SFR. 
The uncertainties are within a factor of about two (e.g., 
Inoue et al. 2000). The other uncertainty is related to 
stellar properties. The choices of the stellar properties 
adopted affect the conversion formula of the SFR. Ken- 
nicutt, Tamblyn, & Congdon (1994) derived the formula, 
SFR/Mq yr-1 = 3.05 x lO-^Lna/i©, adopting their stel- 
lar population synthesis model and Salpeter IMF (0.1- 
IOOM0). We also derive the same type formula by using 
the same IMF but by using very simplified stellar proper- 
ties (Inoue et al. 2001). In comparison with them, we find 
that the uncertainty due to the adopted stellar properties 
is at most a factor of two. Therefore, we conclude that the 
effect of LCE by dust may be more significant than the un- 
certainties due to the specific IMF or the stellar properties 
for nearby spiral galaxies. 

Finally, we discuss the alternative way for determining 
galaxy-wide /. We can use equation (13) to do it. How- 
ever, there are three difficulties to determine / reasonably, 
(i) One is the uncertainty of the cirrus fraction of the IR 
luminosity of the galaxies. It is widely believed that the IR 
luminosity of the galaxies consists of two components: the 
warm component originating from the dust nearby the star 
forming regions and the cirrus component originating from 
the diffuse cold dust heated by the interstellar radiation 
field from the older stars (e.g., Lonsdale-Persson & Helou 
1987). It is a very difficult task to subtract the cirrus com- 
ponent from the observed IR luminosity of the galaxies in 
order to obtain the "pure" IR luminosity originating from 
the star forming regions.*' Since the cirrus fraction differs 
from galaxy to galaxy (e.g., 0.3-0.8; Lonsdale-Persson & 
Helou 1987), the subtraction of cirrus causes a significant 
uncertainty in the estimation of /. 

The other difficulties are the correction of the Ha lu- 
minosity for (ii) the interstellar extinction and (iii) the [N 
11] contamination. Indeed, Kennicutt (1983) noted that 
the interstellar extinction at the wavelength of Ha is 1.1 
mag on average for his sample galaxies. That is, the incre- 
ment correction factor for Ha luminosity is 2.8. For that 
case, we do not neglect such large extinction in order to 
determine / reasonably. At that time, we may also have 
to correct the inclination of the galaxies along the line of 
sight to find the correct /. The [N 11] contamination is 
more complicated if we do not have sufficient resolution of 
the spectrograph. Fortunately, the Ha flux of Caplan et 
al. (2000) is not contaminated by [N 11] because of their 



8 



Inoue, Hirashita and Kamaya 



high dispersion of spectrograph. Then, our main conclu- 
sions on / of the Galactic H ii regions are not altered. 
However, we are sure of that the [N ii] line contaminates 
the galactic Ha emission (e.g. Kennicutt 1983). When we 

use equation (13) to find /, therefore, we always need the 
Ha luminosity corrected for the interstellar reddening and 
[N ii] contamination to find the true /. Unfortunately, we 
face the deficiency of the data to correct the Ha luminosity 
of each individual galaxy. Because of several large uncer- 
tainties, thus, we cannot safely determine / of the galaxies 
via equation (13) in this paper. This problem is examined 
by us in the near future. 

5. CONCLUSIONS 

For each individual H II region in the Galaxy, we exam- 
ine LCE by dust within such regions. Then, we discuss 
the LCE effect on determining the SFR of galaxies. We 
reach the following conclusions. 

[1] We estimate the fraction of Lyman continuum pho- 
tons contributing to hydrogen ionization within the Galac- 
tic H II regions, /. The estimated / values distribute over a 
wide range (0.02-1). In our framework, an H II region with 
higher metallicity and larger Lyman photon flux tends to 
have a smaller /. 

[2] Prom the correlation tests in Figure 3, we find that 
/ is approximately treated as a function of only the dust- 
to-gas mass ratio, V. For our sample of the Galactic H ii 
regions, the Orion nebula is regarded as a typical region 
in gas density and Lyman continuum flux. 



[3] We explain theoretically the observed correlation be- 
tween IR luminosity and radio luminosity of the Galactic 
H II regions. Then, we manifest that a part of Lyman 
continuum photons (about a half or more) ought to be 
absorbed by dust directly in many regions. 

[4] We do not flnd the observational correlation between 
fluxes of IRAS PSC and Ha for our sample H ii regions, 
which is expected by our model theoretically. This may in- 
dicate that IRAS PSC is not suitable to determine the IR 
flux of nearby extended sources like Galactic H II regions 
as pointed out in literatures. 

[5] If we assume that D in H II regions is not signifi- 
cantly different from a mean V of the ISM, we estimate / 
averaged over galactic scale to be about 0.3 for the nearby 
spiral galaxies (i.e. / ~ 0.3). Consequently, the increment 
factor for the SFR due to the LCE by dust (1//) is about 
3. The correction factor is larger than the uncertainties of 
determining the SFR due to the choice of a specific IMF 
or stellar properties (within a factor of 2). Therefore, we 
cannot neglect the effect of LCE by dust when we estimate 
the SFR in galaxies. 
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Table 1 



Representative Giant Galactic H ii Regions 



Object 


Distance 




no 


12 + log (0/H) 


V 


/ 


Reference 




(kpc) 


(1048s-l) 


(cni^'^) 




(10-3) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Orion 


0.5 


7.59 


150 


8.76 


5.0 


0.26 


1 


M8 


1.1 


10.6 


60 


8.74 


4.5 


0.37 


1 


M16 


2.5 


29.1 


(100) 


8.76 


5.0 


0.16 


1 


Rosette 


1.4 


43.2 


13 


8.20 


0.79 


0.85 


2 


M17 


1.8 


143 


(100) 


8.81 


6.3 


0.016 


1 


Carina 


1.4 


167 


13 


8.49 


1.8 


0.56 


2 


NGC3603 


8.5 


897 


25 


8.51 


2.0 


0.23 


2 



Note. — Col. (3): Number of Lyman continuum photons estimated from radio observations by 
using equation (2) in Mezger, Smith, & Churchwell (1974). Col. (4): Averaged electron number 
density from Kennicutt (1984). Two parentheses are assumed values. Col. (5): Observed Oxygen 
abundance from Kennicutt et al. (2000) and references therein. Col. (6): Determined dust-to- 
gas mass ratio by using the ©-(O/H) relation of figure 2 (see also Hirashita 1999a,b). Col. (7): 
Estimated / parameter by using equations (7), and (10). 



References. — Sources of radio observations: 1. Mezger & Henderson (1967); 2. Shaver & Goss 
(1970) 
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Table 2 

Galactic H ii Regions of Caplan et al. (2000) 



Obiect 


-'-'Ha 

{L\j erg s j 


Nl 
(iU s j 


(cm j 


12 + loe- CO/H'l 

\^) 


V 
(iU j 


f 

v) 


S54 


42.0 


3.08 


245 


8.60 


2.7 


0.54 


S88 


2.17 


0.16 


28 


8.43 


1.5 


0.94 


S93 


8.65 


0.63 


269 


8.57 


2.5 


0.71 


SlOl 


3.61 


0.26 


7 


8.34 


1.1 


0.97 


S104 


76.7 


5.62 


56 


8.51 


2.0 


0.71 


S127 


38.4 


2.81 


545 


8.23 


0.79 


0.80 


S131 


0.20 


0.01 


18 


8.39 


1.3 


0.98 


S142 


16.6 


1.22 


21 


8.27 


0.89 


0.94 


S148 


58.2 


4.27 


235 


8.30 


1.0 


0.78 


S153 


8.90 


0.65 


106 


8.12 


0.59 


0.94 


S156 


59.4 


4.35 


907 


8.42 


1.5 


0.55 


S168 


18.4 


1.35 


137 


8.30 


1.0 


0.87 


S184 


4.56 


0.33 


90 


8.39 


1.4 


0.90 


S206 


31.8 


2.33 


313 


8.37 


1.3 


0.75 


S209 


201 


14.7 


645 


8.18 


0.71 


0.69 


S211 


26.0 


1.91 


135 


8.03 


0.48 


0.93 


S219 


4.26 


0.31 


166 


8.27 


0.89 


0.92 


S252 


5.22 


0.38 


144 


8.35 


1.2 


0.89 


S257 


2.55 


0.19 


160 


8.17 


0.69 


0.95 



Note. Col. (2): Ha luminosities are calculated by Lna = 47rD^fHQ, where D is the 
kinematic distance and i^Ha is the Ha flux corrected for interstellar extinction by Balmer 
decrement from Caplan et al. (2000) and the Galactic extinction properties (Savage & 
Mathis 1979). Col. (3): Number of Lyman continuum photons estimated from col.(2). 
Col. (4): Electron number density from Deharveng et al. (2000). Col. (5): Observed 
Oxygen abundance from Deharveng et al. (2000). Col. (6): Dust-to-gas mass ratio 
determined from col. (5) by using the I?-(0/H) relation of figure 2 (see also Hirashita 
1999a,b). Col. (7): Estimated / parameter by using equations (7), and (10). 
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Table 3 



Ha - IR CORRELATION OF THE GALACTIC H II REGIONS 





H/yJD — V J 






FiRAS 


6 


J 




(mag) 


(10-1° erg cm-2) 




(10"* erg 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


S83 


1.75 


6.88 


19223+2041 


0.98 


1 


1.35 


S93 


1.76 


12.6 


19529+2704 


6.91 


1 


0.49 


SlOO 


1.31 


19.4 


19598+3324 


50.8 


1 


0.11 


SlOl 


0.56 


6.07 


19579+3509 


0.22 


0.98 


2.49 


S127 


1.65 


2.77 


21270+5423 


1.29 


1 


0.56 


S128 


1.83 


6.31 


21306+5540 


3.53 


1 


0.48 


S138 


1.38 


2.80 


22308+5812 


4.15 


1 


0.20 


S146 


1.57 


6.16 


22475+5939 


4.05 


1 


0.42 


8148 


0.99 


4.09 


22542+5815 


2.73 


1 


0.41 


S152 


1.10 


9.28 


22566+5830 


6.68 


1 


0.38 


S153 


0.68 


2.45 


22571+5828 


2.84 


0.99 


0.25 


S156 


1.12 


16.5 


23030+5958 


7.34 


1 


0.58 


S168 


0.99 


8.95 


23504+6012 


1.01 


1 


1.55 


S206 


1.36 


52.1 


03595+5110 


2.44 


1 


2.32 


S211 


1.61 


5.08 


04329+5045 


1.29 


1 


0.91 


S212 


0.92 


13.8 


04368+5021 


0.73 


1 


2.22 


S217 


0.73 


2.72 


04551+4755 


0.50 


1 


1.14 


S252 


0.33 


7.90 


06066+2029 


3.82 


0.90 


0.51 



Note. — Col. (2): Color excess determined from the Balmer decrement. Col. (3): Ha fluxes corrected 
for interstellar extinction of col. (2) by the Galactic extinction properties (Savage & Mathis 1979). 
Col.(5): IRAS 40 120 fim fluxes derived by Fjras = 1-26 x 10-"(2-56fco/Jy + -Fioo/Jy)- Col.(6): 
Estimated e by using equation (14). Col.(7): Estimated / by using equation (13) and col. (6). 
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Table 4 

/ OF Seven Nearby Spiral Galaxies 



Object 


Distance 




-P200 


Td 


logMd 


l0gMHI+H2 


V 


/ 




Mpc 


(Jy) 


(Jy) 


(K) 


(Mo) 


(Afo) 


(10-3) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 134 


21.1 


67.0 


127 


21 


8.18 


10.40 


5.9 


0.21 


NGC 628 


8.8 


67.5 


209 


19 


7.72 


10.04 


4.8 


0.29 


NGC 660 


11.3 


120 


158 


24 


7.51 


10.04 


2.9 


0.48 


NGC 5194 


6.2 


303 


533 


22 


7.62 


9.88 


5.5 


0.24 


NGC 5236 


6.9 


624 


889 


23 


7.90 


10.54 


2.3 


0.56 


NGC 6946 


10.1 


338 


743 


20 


8.38 


10.48 


8.0 


0.12 


NGC 7331 


10.9 


120 


243 


21 


7.85 


10.00 


7.0 


0.16 



Note. — Col.(l): Object name. Col. (2): Distance from our Galaxy in Mpc from Alton 

et al. (1998). Cols.(3) - (4): Integrated flux density of ffiylS-HIRES image at 100, 200 
liui, respectively, from Alton et al. (1998). Col. (5): Dust temperature estimated from 
F1Q0/F200 and the modified black-body spectrum with the spectral index of 1. Col. (6): 
Dust mass estimated from i^ioo and Ta via equation (15). Col. (7): Total (H i + H2) 
gas mass. The data of six galaxies except NGC 134 are taken from Devereux & Young 
(1990). H I and H2 data of NGC 134 are taken from Mathcwson & Ford (1996) and 
Elfhag et al. (1996), respectively. In the calculation of mass, the same CO-H2 conversion 
factor, Xco = 2.8 x 10^° [(H2 cm-^)/(K km s"^)], and the distances in col.(2) are used 
for all the sample galaxies. Col. (8): Dust-to-gas mass ratio. Col. (9): / estimated V and 
equations (10) and (11) with x = 2. 




Fig. 1. — / vs. dust optical depth over actual ionized radius, ra- 
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Fig. 2. — Global relation between V and Oxygen abundance proposed by Hirashita (1999a,b). 
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Fig. 3. — Various correlations among various properties of our sample of the Galactic H ii regions. Open symbols are the samples in Table 1, 
and Filled symbols are the samples in Table 2. (a) metallicity vs. observed Lyman continuum number per unit time, N^^. (b) metallicity vs. 
electron number density, Tig. (c) Ni vs. Ue- (d) Estimated / vs. metallicity. (e) / vs. JV^y- (f) / vs. ng. 




Fig. 4. — Estimated / vs. dust-to-gas mass ratio, V. Open symbols are the representative giant H II regions in Table 1, and filled symbols 
are H ii regions in Table 2. This figure is basically the same as figure 3 (d). The solid line denotes the best fit, x = 1.8 in equation (11). Two 
dotted lines mean x = 1.0 and 3.0. 
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Fig. 5. IRAS flux vs. Hq flux of some Galactic H II regions. The plotted points are the sample H II regions of Table 3. The lines are the 
Piras^Fhci relation expected by equation (13). The solid, dotted, and dashed lines mean / = 1, 0.5, and 0.1, respectively, adopting e = 1. 
We set / = 1, e = for the dash-dotted line. 
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log IR luminosity (solar unit) 




Fig. 6. — Relation between IR luminosity and number of Lyman continuum photons. The crosses are observational data for H II regions 
reported by Wynn- Williams & Becklin (1974). The parameter sets of lines are the same as those of figure 5. 



